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We find that recent results from the KamLAND collaboration on geologically produced antineu- 



trinos, N(U+Th) = 28l]J events, correspond to a radiogenic heat production from Uranium and 
& Thorium decay chains H(U+Th) = 38^33 TW. The 99% confidence limit on the geo-neutrino signal 



translates into the upper bound H(U+Th) < 162 TW, which is much weaker than that claimed 



by KamLAND, H(U+Th) < 60 TW, based on a too narrow class of geological models. We also 
performed an analysis of KamLAND data including recent high precision measurements of the 
13 C(a,n) 16 cross section. The result, N(U+Th) = 3ll} 3 , corroborates the evidence (~ 2.5a) for 
geo-neutrinos in KamLAND data. 
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I. INTRODUCTION 

Geologically produced antineutrinos (geo-neutrinos) were introduced by Eder [l| in the 
sixties and Marx soon realized their relevance. In the eighties Krauss et al. discussed their 
potential as probes of the Earth's interior in an extensive publication ^|. In the nineties 
the first paper on a geophysical journal was published by Kobayashi et al. Q]. In 1998, 
Raghavan et al. \^ and Rotschild et al. [f| pointed out that KamLAND and Borexino 
should be capable of geo-neutrino detection. 

The potential of geo-neutrinos for providing information on the energetics and compo- 
sition of the Earth has been discussed in Refs. [7, 0] and more recently in Ref. The 
KamLAND collaboration has just published Q| their experimental results, claiming some 
28 geo-neutrino events from Uranium and Thorium decay chains in a two-year exposure. 
This important step shows that the technique for exploiting geo-neutrinos in the investi- 
gation of the Earth's interior is now available. In order to understand where to go with 
geo-neutrinos, we have to know where we stand, in the light of the available data. In this 
spirit, the aim of this letter is to discuss the implication of the KamLAND result on the 
contribution of Uranium and Thorium decay chains to the terrestrial heat. 



II. THE GEO-NEUTRINO SIGNAL AND THE RADIOGENIC TERRESTRIAL 
HEAT 

For a given value of, e.g., the Uranium mass in the Earth, m(U), the contributed heat 
production rate from the Uranium decay chain is uniquely determined, H(U) = 0.95 TW 
xm(U)/(10 16 kg), whereas the flux and signal of geo-neutrinos depend on the detector loca- 
tion and on the Uranium distribution inside the Earth. The connection between the signal 
of geo-neutrinos from the Uranium decay chain, the mass of Uranium in the Earth and the 
heat production rate from that element was found in Ref. |9J, by using global mass bal- 
ance together with a detailed geochemical and geophysical study of the region surrounding 
Kamioka. 

We remark that the mass of Uranium in the crust, m c (U), is rather well constrained by 
geological data, in the interval (3-f 4) x 10 16 kg. The main uncertainty is the amount of 
Uranium (and Thorium) in the mantle. Geo-neutrinos should provide us with this informa- 
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tion. 

For a total uranium mass m(U) in the Earth, the maximal and minimal signals can be 
derived by using a proximity argument: the maximal (minimal) signal is obtained by placing 
the sources as close (as far) as possible to (from) the detector, consistently with geochemical 
and geophysical constraints. The maximal signal is thus obtained for an Uranium rich 
crust, m c (U) = 4 x 10 16 kg, and distributing uniformly in the mantle the rest of Uranium, 
m(U) — m c (U). The minimal signal corresponds to an Uranium poor crust, m c (U) = 3 x 10 16 



re lower mantle. These maximal and 
where the interested reader can find 



OllC 



kg, and distributing the rest near the bottom of t 
minimal signals provide the borders of Fig. 5 of Ref 
more details. 

We have extended the same calculations to Thorium, assuming a global chondritic 
Uranium-to-Thorium mass ratio, m(Th)/m(U) = 3.9 ± 0.1, so that we can now connect 
the combined signal of geo- neutrinos from Uranium and Thorium progenies, S(U+Th), with 
the radiogenic heat production rate from these elements, H(U+Th), see Fig. ^ 

The geo-neutrino signal is expressed in Terrestrial Neutrino Units, one TNU correspond- 
ing to 10~ 32 z7 e captures per target proton per year. 

The allowed band in Fig. ^ is estimated by considering rather extreme 1 models for the 
distributions of radioactive elements, chosen so as to maximize or minimize the signal for a 
given heat productions 3. 

We also remark that, in comparison with the experimental error, the width of the band is 
so narrow that we can limit the discussion to the median line of the allowed band in Fig. CJ 
which represents our best estimate for the relationship between signal and power. 

For the Bulk Silicate Earth (BSE) model, H(U+Th) = 16 TW, our prediction for Kamioka 
is centered at 37 TNU. 

By assuming that Uranium and Potassium in the Earth are in the ratio 1/10,000 and 
that there is no Potassium in the core, the total radiogenic power is H(U+Th+K) = 1.18 
H(U+Th). With these assumptions, a maximal and fully radiogenic heat production rate, 
H(U+Th+K) = 44 TW, corresponds to H(U+Th) = 37 TW, which gives a signal S(U+Th) 
w 56 TNU. 



1 We are neglecting here the possibility that some Uranium or Thorium is hidden in the core. This possibility, 
which has been advanced by some authors (see, e.g., Herndon 11] and Hofmeister and Criss [l^) would 
imply an even smaller signal-to-power ratio. 
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The KamLAND collaboration has reported Q| data from an exposure of N p = (0.346 ± 
0.017) x 10 32 free protons over a time T = 749 days with a detection efficiency e = 69%: 
the effective exposure is thus E eff =N p x T x e = (0.487 ± 0.02 5) x 10 32 protons • yr. In 
the energy region where geo-neutrinos are expected, there are C = 152 counts, implying a 
statistical fluctuation of ±12.5. Of these counts, a number R = 80.4 ± 7.2 are attributed to 
reactor events, based on an independent analysis of higher energy data. Fake geo-neutrino 
events, originating from 13 C(a, n) w O reactions following the alpha decay of contaminant 
210 Po, are estimated to be F = 42 ± 11, where the error is due to a 20% uncertainty on 
the 13 C(a,n) 16 cross section and a 14% uncertainty on the number of 210 Po decays in the 
detector. Other minor backgrounds account for B = 4.6 ± 0.2 events. The number of geo- 
neutrino events is estimated by subtraction, N(U+Th) = C- R- F-B, with an uncertainty 
obtained by combining the independent errors: N(U+Th) = 25ljg. The geo-neutrino signal 
is thus S(U+Th) = N(U+Th) / E cfr = 51+j| TNU. From the median line in Fig. done finds 

H(U+Th) = 3llg TW (rate only) . (1) 

This "rate only" study has been improved in Ref. ^(| by exploiting the shape of the spec- 
trum. A likelihood analysis of the unbinned spectrum yields N(U+Th) = 28l]J, see Fig. 4b 
of Ref. Q. This implies S(U+Th) = 57±|f TNU and 

H(U+Th) = 3811 TW ( rate + spectrum) . (2) 

The best fit value is close to the maximal and fully radiogenic model, however the BSE is 
within la. 

By using the median line in Fig. ^ the 99% confidence limit on the signal (145 TNU) 
corresponds to 133 TW. If we include the uncertainty band of the theoretical models, we 
find an upper bound of 162 TW, see point A in Fig. ^ This point corresponds to a model 
with a total Uranim mass ra(U) = 80 x 10 16 kg, an uranium poor crust, m c (U) = 3 x 10 16 
kg, the rest of the Uranium being placed at the bottom of the mantle, and global chondritic 
Thorium-to-Uranium ratio. 

This 162 TW upper bound is much higher than the 60 TW upper bound claimed in 



Ref. 
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, which was obtained by using a family of geological models which are too narrow 



and are also incompatible with well-known geochemical and geophysical data. 
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In fact, the authors of Ref. [10] start with a reference BSE model derived from Ref. jl^ : 
the total Uranium mass is mesE = 8 x 10 16 kg, roughly half in the crust and the rest in the 
mantle, and the abundance ratio is Th/U = 3.9. This model corresponds to Ubse(U + T1i) = 
16 TW and predicts a signal of 38.5 TNU, very close to our prediction for BSE. The signal 
is assumed to scale with the total mass of U+Th, so that heat production and signal are 
also proportional: 

S(U+Th) = 38.5 TNU x H(U+Th)/(16 TW) . (3) 

In this way, the 99% upper limit on the signal, 145 TNU, is translated into 60 TW 4(J, see 
point B in Fig. ^ 

This scaling assumption, however, produces a too limited series of models. The points 
in the shaded area of Fig. ^ correspond to all models which are compatible with available 
geochemical and geophysical data 2 : most of these models cannot be obtained by Eq. Eland 
predict, for a given signal, a larger power than Eq. El which therefore cannot be used to 
derive an upper bound on the radiogenic power production. 

Furthermore, Eq. El implies that Uranium in the crust (and in the mantle) scales linearly 
with the total Uranium mass. This becomes incompatible with the geochemical data on the 
crust (m c (U) < 4 x 10 16 kg) already for total masses slightly above the BSE estimate, i.e., 
for models where H(U+Th)> 20 TW. For example, in the model yielding 60 TW (point B 
in Fig. the crust should contains about 13 x 10 16 kg of Uranium, four times more than 
the largest geochemical estimate. This inconsistency is clearly seen in Fig. which shows 
that the family of models labeled as "rescaled models" lies essentially in the geo-chemically 
excluded region. 

We remark that the bound H(U+Th)<162 TW which we have extracted from KamLAND 
data does not add any significant information on Earth's interior, since anything exceeding 
H(U+Th)= 37 TW (i.e. H(U+Th+K)=44 TW) is unrealistic. The upper limit simply 
reflects the large uncertainty in this pioneering experiment. 

On the other hand, what is important for deciding the potential of future experiments is 
the relationship between geo-neutrino signal and heat production in the physically interesting 
region, H(U+Th)< 37 TW. The basic parameter is the slope, dS/dH, which expresses how 



2 We note that actually models with H(U+Th)> 37 TW are essentially unrealistic; this portion of the graph 
is included just for discussing KamLAND results. 
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the experimental error translates into an uncertainty on the deduced heat production. For 
our models we find from Fig. ^dS/dH ~ 1 TNU/TW. Discrimination between BSE and 
fully radiogenic models, which requires a precision AH ~ 7 TW, requires thus an experiment 
with an accuracy AS ~ 7 TNU. The "rescaled models" of ref. [10j, yielding dS/dH~ 2.4 
TNU/TW, might misleadingly suggest that the same goal can be achieved for AS=17 TNU. 



uncertainty has been adopted in 
experimental papers (see, e.g., Ref. 



III. THE GEO-NEUTRINO SIGNAL AND THE 13 C(a,ra) 16 CROSS SECTION. 

As already remarked, a major uncertainty for extracting the geo-neutrino signal originates 
rom the 13 C(a, n) 16 cross section 3 . The values used in Ref. [lG| are taken from the JENDL 
1^ | compilation, which provides an R-matrix fit of relatively old data. A 20% overall 

, corresponding to the accuracy claimed in the original 

m 

Recently a series of high precision measurements for this cross section has been performed 
In the relevant energy range (1 -j- 5.3) MeV, the absolute normalization has been 
determined within a 4% accuracy. The measured values are generally in very good agreement 
with those recommended in JENDL, see Fig. El however, we find that the neutron yield per 
alpha particle is 5% smaller. It follows that the number of fake neutrinos is lower, F = 
40 ± 5.8, and geo- neutrino events obviously increase 4 . 

The "rate only" analysis gives now 271^ geo-neutrino events, corresponding to S(U+Th) 
= 551^ TNU. From the median line of Fig. d the radiogenic power is now: 

H(U+Th) = 36tH TW (rate + new 13 C(a,n) 16 0) . (4) 

We also performed an analysis 5 of the binned spectrum reported in Fig. 3 of Ref. [lo| . 
This analysis gives N(U+Th) = 3lt\t counts, corresponding to S(U+Th) = 63±| TNU and 



3 In fact, the claim of 9 geo-neutrino events in Ref. should be dismissed: more than half of these events 
are to be considered as fake signal, produced from 13 C(a, n) 16 reaction. 

4 Indeed Ref. jlCJ mentions that an alternative analysis including the time structure of the scintillation light 
from different particles produced a slightly larger geo-neutrino signal, which is consistent with the result 
presented here. 

5 A complete analysis requires several details (the un-binned spectrum, the energy dependence of the de- 
tection efficiency, . . . ) which are not available to us. Just for a comparison, the binned spectrum analysis 
using the JENDL cross sections with 20% uncertainty gives us N(U+Th)=28.5^4, in agreement with 
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thus: 

H(U+Th) = 44+^ TW (rate + spectrum + new 13 C(a,n) 16 0) . (5) 

IV. CONCLUDING REMARKS 

In summary, the new data on 13 C(a,n) 16 corroborate the evidence for geo-neutrinos in 
KamLAND data, which becomes near to 2.5a. 

On the other hand, the determination of radiogenic heat power from geo-neutrino mea- 
surements is still affected by a 70% uncertainty. The best fit of H(U+Th) is close to the 
prediction of a maximal and fully radiogenic model, however the BSE prediction is within 
la from it. 

With more statistics KamLAND should be capable of providing a three sigma evidence 
of geo-neutrinos, but discrimination between BSE and fully radiogenic models definitely 
requires new detectors, with class and size similar to that of KamLAND, far away from 
nuclear power plants. 

Acknowledgments 

We are grateful to C. Rolfs and his group for useful discussions and for allowing us to 
use their results. 

We thank for their useful comments A. Bottino, E. Lisi, W. F. McDonough, and R. Ragha- 
van. We appreciated the suggestions of the anonymous referee. 



[1] G. Eder, Nucl. Phys. 78 (1966) 657. 

[2] G. Marx, Czech. J. Phys B 19 (1969) 1471. 

[3] L.M. Krauss et al., Nature 310 (1984) 191. 

[4] M. Kobayashi and Y. Fukao, Geophy. Res. Lett. 18 (1991) 633. 

[5] R.S. Raghavan et al., Phys. Rev. Lett. 80 (1998) 635. 

[6] GC. Rotschild et al., Geophy. Res. Lett. 25 (1998) 1083. 

[7] G. Fiorentini, F. Mantovani and B. Ricci, Phys. Lett. B 557 (2003) 139 
|arXiv:nucl-ex/0212008l . 



7 



[8] F. Mantovani, L. Carmignani, G. Fiorentini and M. Lissia, Phys. Rev. D 69 (2004) 013001 
|arXiv:hep-ph/03090l"3] . 

[9] G. Fiorentini, M. Lissia, F. Mantovani and R. Vannucci, Phys. Rev. D 72 (2005) 033017 

|arXiv:hep-ph/0501lTT] . 
[10] T. Araki et al (KamLAND coll.), Nature 436 (2005) 499. 
[11] J. M. Herndon, Proc. Natl. Acad. Sci. USA 93 (1996) 646-648, Geophysics. 
[12] A.M. Hofmeister and R.E. Criss, Tectonophysics 395 (2005) 195. 
[13] S. Enomoto, Thesis, Tohoku Univ. (2005), 

http : / / www. awa.tohoku.j p /KamLAND / publications /Sanshiro_thesis .pdf , 
[14] JENDL Japanese Evaluated Nuclear Data Library, 



http:/ /wwwndc. tokai.jaeri.go.jp/jendl/ 



[15] S. Harissopulos et al., Phys. Rev. C (2005) to appear arXiv:nucl-ex/0509014 . 

[16] K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90 (2003) 021802 

|arXiv:hep-ex/021202l"1 . 
[17] J.K. Bair and F.X. Haas Phys. Rev. C 7 (1973) 1356. 



8 




FIG. 1: The combined signal from Uranium and Thorium geo- neutrinos and the radiogenic heat 
production rate. 

The shaded area denotes the region allowed by geochemical and geophysical constraints. The 
dashed median line represents our best estimate for the relationship between signal and power. 

The dotted line denotes the "rescaled models" of Eq. |3J used in 10!]. Note that most of these 
models are outside the allowed area. 

One TNU corresponds to 10~ 32 v e captures per target proton per year. 
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FIG. 2: Cross section of 13 C(a,n) 16 0. The solid line corresponds to the JENDL compilation, dots 
are the experimental points from Ref. 
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